Adiabatic transfer of light in a double cavity and the optical Landau-Zener problem 
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We analyze the evolution of an electromagnetic field inside a double cavity when the difference 
in length between the two cavities is changed, e.g. by translating the common mirror. We find 
that this allows photons to be moved deterministically from one cavity to the other. We are able 
to obtain the conditions for adiabatic transfer by first mapping the Maxwell wave equation for the 
electric field onto a Schrodinger-like wave equation, and then using the Landau-Zener result for the 
transition probability at an avoided crossing. Our analysis reveals that this mapping only rigorously 
holds when the two cavities are weakly coupled (i.e. in the regime of a highly reflective common 
mirror), and that, generally speaking, care is required when attempting a hamiltonian description 
of cavity electrodynamics with time-dependent boundary conditions. 



I. INTRODUCTION 

One of the outstanding challenges facing quantum in- 
formation proposals based on cavity quantum electro- 
dynamics (CQED) is the development of a method for 
transferring photons between different cavities [T]. With 
a robust photon transfer method, a cavity containing one 
or more atoms can act as a single node in a quantum net- 
work made up of many such cavities connected by optical 
fibres Quantum networks provide a way to solve the 
problem of scalability in quantum information process- 
ing, scalability being one of the Divincenzo criteria for 
the realization of a quantum computer [3] . 

A scheme for the efficient transfer of photons between 
distant atoms has been put forward by Cirac et al. [3]. 
They considered what happens when a photon wave 
packet is emitted by an atom into a high-Q optical cav- 
ity which is connected, via an optical fibre, to a second 
cavity containing a second atom. A generic wave packet 
reflects from the highly reflective mirror of the second 
cavity, thus leaving the wave packet to slosh around in- 
side the cavity-fibre system. However, Cirac et al. real- 
ized that if the emitted wave packet was time-symmetric 
then it could be absorbed by the second atom with high 
efficiency. The coherent exchange of photons between 
an atom and a cavity mode can be controlled by time- 
dependent laser pulses applied through the side of a cav- 
ity [S], and this allows the temporal shape of the wave 
packet to be controlled. 

A different approach has been suggested by Pellizari 
[6l |7j . In analogy with stimulated Raman adiabatic pas- 
sage (STIRAP) [51 H], which uses two or more control 
lasers to transfer an atom between two internal states 
via an intermediate state which is never populated, Pel- 
lizari proposed that laser pulses applied to the two atoms 
could be used to control the transfer of photons between 
them along the fibre, the role of the intermediate state 
being played by the cavity and fibre photon states. The 
appeal of an adiabatic passage method lies in the fact 
that the details of the pulse shape and timing are not 
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Figure 1. (Color online) Double cavity setup consisting of two 
perfectly reflecting mirrors separated by a partially transmis- 
sive common central mirror. AL = L\ — L2 is the difference 
in length between the two cavities. 



important. Furthermore, because the cavity and fibre 
modes are only virtually populated, the scheme is insen- 
sitive to cavity and fibre losses. Inspired by [6], the theory 
of adiabatic passage of photons between two atoms via 
cavities and optical fibres has been developed by a num- 
ber of groups, see e.g., [im - [T2] . We also note that the 
basic concept of STIRAP has also been discussed in the 
context of the adiabatic passage of matter waves, such 
as atoms and electrons, where it is referred to as coher- 
ent tunnehng adiabatic passage (CTAP). For example, 
references [T3] and [H] consider the adiabatic transfer of 
ultracold atoms along a chain of traps. 

In this paper we consider the adiabatic transfer of light 
between two coupled cavities by changing their lengths. 
This can be achieved, for example, by translating the cav- 
ity mirrors using piezo-electric actuators (motors). Al- 
ternatively, if the cavity takes the form of a dielectric 
waveguide the mirror positions can remain fixed but the 
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optical length of the cavities can be controlled by mod- 
ulating the refractive index [15]. The specific model sys- 
tem we choose to illustrate the adiabatic transfer process 
in either of these experimental scenarios is a double cav- 
ity in one dimension as depicted in Figure [l] Changing 
the lengths of the two cavities changes the structure of 
the electromagnetic modes in such a way that a mode 
which is initially localized in one cavity can be smoothly 
transferred into the other cavity. Our study of the spa- 
tial adiabatic passage of light is in the same spirit as 
Pelizzari's, although we shall not consider the atoms, but 
focus instead on the global normal modes of Maxwell's 
equations whose frequencies form a net of avoided cross- 
ings as a function of mirror displacement. It would be 
straightforward to introduce atoms into the scheme: for 
certain mirror displacements a quasi-resonant global light 
mode can be strongly localized in one of the cavities, and 
so a photon emitted by an atom in that cavity will find 
itself in that chosen mode. This global mode is then 
transferred to the other cavity, where the other atom 
can interact with the photon with high efficiency. On 
the other hand, if during the emission the mirror is not 
at a position which localizes the relevant global mode 
in the first cavity, then the photon will be placed into 
a superposition of global modes and the resulting time- 
dependence of the superposition leads to an oscillation 
of the photon back and forth between the two cavities, 
thereby lowering the transfer efficiency to the other atom 
and/or introducing the necessity for good timing. 

Because our model is simple, it allows a straightfor- 
ward analysis of the conditions for adiabatic transfer, 
and a number of analytic results can be obtained. Fur- 
thermore, the Yale opto-mechanical group have realized 
a double cavity (TBI [IT] which is close to that consid- 
ered by us here. Their experimental setup consists of 
an optical cavity divided in half by a thin silicon mem- 
brane. When the cavity is driven through an end mirror 
by an external laser, the membrane vibrates due to ra- 
diation pressure. Indeed, a recent theoretical paper |18j . 
discusses photon shuttling between the two halves of the 
cavity, in close correspondence to what we will consider 
here. Other theoretical studies of this system treat it as a 
photonic version of the Josephson junction ir^HH]. The 
important difference between our paper and previous op- 
tomechanical studies is that we propose to control the 
motion of mirror from outside rather than allowing it to 
move under the action of radiation pressure. 

The problem of a cavity with moving mirrors is also of 
general interest from the point of view of studying CQED 
with time-dependent boundary conditions. This is an 
interesting generalization of CQED with many poten- 
tial applications, including the field of opto-mechanics. 
Other related applications include the possibility of mod- 
ifying the decay rate of light from a cavity by modulating 
the position of the mirrors in time [22H24] , and if the mir- 
rors accelerate very rapidly, photon pair production [25l - 
[27] . Observations of this dynamical Casimir effect have 
recently been reported in microwave cavities |28| . Also 



on the experimental side, advances in the ultrafast ma- 
nipulation of the refractive index of silicon microcavities 
(which is equivalent to changing the cavity's length) have 
been used to demonstrate transitions between modes ^5S] . 
When the refractive index is changed more slowly, the 
light adiabatically follows the modes, leading to a change 
in its color [Tl|[3nHS2- 

The programme we follow in this paper is to treat cav- 
ity electrodynamics with time-dependent boundary con- 
ditions from first principles, by starting with Maxwell's 
wave equation. From this we derive an effective 'hamil- 
tonian' which can be used for the purposes of first-order 
time propagation. More precisely, by making an approxi- 
mation analogous to the paraxial approximation familiar 
from optics, but in time rather than space, we find that 
the Maxwell wave equation reduces to an equation of 
motion for the classical Maxwell field that is first order 
in time and has a mathematical resemblance to the time- 
dependent Schrodinger equation |22j . (For a different ap- 
proach to mapping the Maxwell wave equation onto the 
Schrodinger equation see [53] [51]). Under certain condi- 
tions we are therefore able to apply Landau-Zener theory 
|35H37j ■ familiar from non-relativistic quantum mechan- 
ics, to the calculation of the coupling between classical 
Maxwell modes induced by changing the lengths of the 
cavities. We emphasize that our method is different from 
one in which the hamiltonian obtained from the static 
problem is used for time propagation. In particular, we 
show that there are higher order terms which provide 
significant corrections to the latter scheme [52]. 

The idea of using the Landau-Zener model to under- 
stand time-dependent processes in optics is not new. For 
example, the experiment reported in |38j studied the time 
evolution of light in cavity where an electo-optic modu- 
lator was used to couple the two polarization states of a 
single longitudinal mode. A second electo-optic modula- 
tor was used to vary the energy separation between two 
polarization states: varying this separation linearly in 
time whilst maintaining a constant coupling corresponds 
to the Landau-Zener model. Another system which is 
useful for studying the analogies between optical and 
quantum dynamics is provided by evanescently coupled 
waveguides, where the coupling between the waveguides 
changes with distance along the waveguide [39]. This 
system has been studied experimentally by a number of 
groups and various optical analogues to quantum phe- 
nomena have been observed, including Zener tunneling 
j41j . adiabatic evolution [40], and dynamical localization 
[32] ■ There have also been numerous studies of the condi- 
tions necessary for adiabaticity in CQED when atoms are 
included, see, for example, [32]. In [33], the dynamics of a 
two-level system coupled to a cavity field is studied in the 
situation where the separation between the two energy 
levels is varied linearly in time, and in |45j the principle 
of STIRAP is applied to a two-mode Jaynes-Cummings 
model with degenerate mode frequencies, with the aim of 
adiabatically transferring photons between modes. We 
show in this paper that the application of the Landau- 
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Zener model to the problem of a double cavity with 
changing cavity lengths illustrates in a simple fashion 
some of the similarities and differences between classical 
electrodynamics obeying Maxwell's equations and quan- 
tum mechanics obeying Schrodinger's equation. 

This paper is organized as follows. In Section [IT] wc 
introduce a simple model for the spatial dependence of 
the dielectric permittivity function inside a double cavity. 
This model treats the central mirror as a Dirac (5-function 
which facilitates analytic calculations later in the paper. 



In Section III we find the global static solutions (normal 



modes) of Maxwell's wave equation subject to this di- 
electric function. In Section |TV] we use these solutions to 
compute the transfer ratio, i.e. the degree to which the 
light can be localized on one side of the central mirror 
or the other depending on mirror position. In Section 
|V]we introduce time-dependence by relating the problem 
of the transfer of light in a double cavity with a mov- 
ing mirror to the well known Landau-Zener problem in 
quantum mechanics. In order to apply the results from 
the Landau-Zener problem we need to map the Maxwell 
wave equation, which is second order in time, onto the 
Schrodinger wave equation, which is first order in time, 
and this is accomplished in Section \VJ\ The electric field 
is quantized in Section jVIlj and in Section [VIII| we ex- 
amine the regimes of validity of the mapping onto the 
Schrodinger wave equation. The (5-function mirror model 
is compared to a more realistic model for a mirror in Sec- 
tion |IX[ and in Section [X] we discuss the experimental 
feasibility. Conclusions are drawn in Section [XI 

This paper has six appendices: in Appendix]^ we com- 
pare and contrast the double cavity problem for light 
with the double well problem in quantum mechanics. In 
Appendix [B] we give some results on the transmission 
amplitudes for finite and ^-function mirrors, and in Ap- 
pendix [C] we explain how the results given in this paper 
can be applied to the case of two coupled waveguides with 
a fixed mirror position but two refractive indices, ni{t) 
and n2{t), that can be separately varied in time. In Ap- 
pcndix[D]we define the diabatic basis and explain its con- 
nection to the adiabatic (global) basis. In Appendix[E|we 
discuss the approximation whereby the time-dependence 
of diabatic basis functions is neglected. In Appendix [F] 
we briefly consider modifications to Maxwell's wave equa- 
tion due to moving dielectric mirrors. 



II. (5-FUNCTION MIRROR MODEL 

Consider a double cavity formed from two end mirrors 
plus a common mirror located between them, as shown 
Figure [l] We assume that the end mirrors are much 
more reflective than the central mirror and that the loss 
of light from the double cavity is negligible during the 
time taken for the light to be transferred from one side 
to the other. The stringent conditions this places on an 
experiment will be examined in Section [Xj As mentioned 
above, our proposed setup is motivated by remarkable 



developments in opto-mechanics [111 HZ] • The additional 
feature we require here is that the common mirror be 
externally controlled, i.e. can be moved along the cavity 
axis. 

A simple theoretical model describing a double cavity 
has been given in a classic paper by Lang, Scully and 
Lamb [15]. For the purposes of solving Maxwell's wave 
equation in the double cavity, they treated the end mir- 
rors as perfect reflectors and the central mirror as a thin 
slab of dielectric material which is modelled by a Dirac 
(5-function spatial profile. The double cavity model is 
thereby encoded in a dielectric permittivity function of 
the form 
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where x — —Li, and x = L2 are the positions of the 
end mirrors, a is a parameter which determines the re- 
fiectivity of the common mirror, see the Appendices for 
more details. The total length of the double cavity is 
L = Li + L2, and we also define the difference between 
the lengths of the two cavities to be AL = L1 — L2, which 
is also twice the displacement of the common mirror from 
the center of the whole cavity. 

Maxwell's wave equation for the electric field £{x, t) in 
the double cavity is 



d^£ix,t) 
dx"^ 



- + aS{x)) 



d^£(x,t) 



= 



(2) 



We use this (5-mirror model for the bulk of this paper 
because its simplicity facilitates analytic results. How- 
ever, in Section IX we compare the results of the 5- 



mirror model to the more realistic case of a mirror of 
finite width. 

Note that for mathematical convenience we have cho- 
sen a coordinate system in Eq. ([T]) where the common 
mirror is always located at x = 0. However, we do not 
intend to attach any physical significance to this choice: 
physically, one can either take the view that the com- 
mon mirror has a fixed position and it is the two end 
mirrors that are displaced, or vice versa. The latter view 
corresponds to an experimental situation, such as that 
described in |16| and |17j , where a common membrane is 
displaced, whilst the former view is closer in spirit to an 
experiment in which the position of the common mirror 
is held fixed, but the refractive index of the two cavities 
is modulated [15], see Appendix [C] for more details. In 
Sections |V] and VI we treat the case of moving mirrors: 
providing the mirrors are moving at constant velocities, it 
is physically equivalent to imagine either translating the 
end mirrors or the common mirror because these two situ- 
ations simply correspond to two different inertial frames. 
However, there are some interesting considerations to be 
taken into account when treating light propagating in 
moving dielectrics, as briefly discussed in Appendix [F] 
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III. SOLUTIONS OF MAXWELL'S WAVE 
EQUATION 

We write the solutions to the Maxwell wave equation as 
£m{x,t) = Umix)eicp{-iujmt), where = fcm/V^oMo is 
the angular frequency and m = 1, 2, 3 . . . is an integer 
labelling the modes. The dimensionless mode functions 
Um{x) can be chosen to be orthonormal in the Sturm- 
Liouville sense by ensuring that they obey 
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£{x)Ul{x)U,n{x)dx = Slr^ 



^0 J-Li 

Inserting the above form for £{x,t) into Eq. ^ gives 

d^Urnix) 



(3) 



+ k'„,il + aSix))Uraix) ^ 



(4) 



Solutions satisfying the boundary conditions Um{—Li) 
Um{L2) = are given by 



Am sin [k 

n 



B„iSin [km{x - Li)\ 



-Li < x < 
< X < La 



(5) 



One boundary condition is given by the continuity of the 
electric field across the (5-mirror, ?7m(0"'") — J7m(0~). A 
second boundary condition is given by integrating Eq. Q 
over a vanishingly small interval containing the mirror, 
leading to C/;„(0+) - U'^^^-) = -a/c^C/™(0). 

Combining the two boundary conditions one is led to 
the following equation for the wave numbers km of the 
allowed modes .46. 



tan(fcmL2) 



tan(A:,„Li) 



afcm tan(fc,„ii) — 1 



(6) 



This transcendental equation can in general only be 
solved numerically. It is useful to re-write it as 



cos{km^L) — cos{kmL) = 2 



sin km L 



(7) 



so that when ak is large the sine function on the right 
hand side is small. The left hand side may then be ex- 
panded around its roots and this permits approximate 
analytic solutions (Eqns ^ and ( |lO| ) below). 

Consider first the situation when the common mirror is 
perfectly refiecting (a — >■ oo), so that the two sides of the 
double cavity are uncoupled. The modes in this case are 
specified by the wave numbers kn = Iwe j {L±AL), which 
come in pairs that cross at AL = 0, and are plotted as 
the dashed red lines in Figures [2] and |3] The mode whose 
wave number decreases with increasing AL corresponds 
to light trapped on the left hand side of the mirror, and 
vice versa. When the modes are coupled there are still 
two for each value of the integer n, whereas each value of 
the integer m introduced previously labels a single mode: 
from henceforth we shall use the "n" convention. 

We refer to the solutions of the wave equation in the 
perfectly reflecting case as localized modes. To find the 
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Figure 2. (Color online) The wave numbers k of the modes 
in a double cavity plotted as a function of AL = L\ — L2 
which is the difference in length between the two cavities. 
The wave numbers are obtained by solving the transcendental 
equation ([7|. When the reflectivity of the common mirror is 
inflnite (a — ^ 00), we obtain the localized modes (dashed 
red lines), which are trapped on one side or the other of the 
mirror. When the reflectivity is finite we obtain the global 
modes (solid blue curves) which feature avoided crossings. 
The total cavity length is set at L = 1 x 10~'*m, and the 
global curves have q = 1 x 10~^m, which translates to a 
transmission probability of the common mirror of about 4%. 



modes when the coupling between the two sides of the 
cavity is finite, one must solve Eq. ([t]) for fc, with a ap- 
propriately chosen to give the desired reflectivity of the 
common mirror. As shown in Appendix [B] the reflection 
probability from the (5-mirror is given by 



R 



k^o?/A 
1 + fc2a2/4 



(8) 



We shall refer to solutions of Eq. ([7| for finite a as global 
solutions, because they extend throughout the double 
cavity, and they are plotted as the solid curves in Fig- 
ures [2] and |3] As can be seen in these figures, the effect 
of the coupling is to turn crossings into avoided crossings. 

At AL = the cavity is symmetric about the central 
mirror and so the solutions of Maxwell's wave equation 
must have well defined parity: for each n we find a pair 
of modes, one of which is even and the other odd. As 
the common mirror is moved away from the center of the 
cavity the global modes lose their well defined parity, but 
we shall continue to label them as "even" and "odd" for 
the sake of simplicity. By expanding Eq. ([7| to second 
order in AL about AL = 0, we obtain the following ex- 
pressions for the wave numbers of the global modes valid 
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and (10), respectively. We find 
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Figure 3. (Color online) A zoom-in of one of the avoided 
crossings shown in Figure [2] corresponding to n = 128. Three 
pairs of curves are shown: the exact numerical solutions of 
Eq. ^ are plotted as solid (blue) curves, the approximate 



analytical solutions given by Eqns (111 and (12 1 are plotted 
as short dashed (black) curves, and the localized solutions are 
plotted as long dashed (red) curves. The difference between 
the exact numerical solutions and the approximate analytical 
ones is hard to discern. Note that in contrast to the double 
well problem in quantum mechanics, the even mode lies at a 
higher frequency uj = ck than the odd mode. 



close to avoided crossings centered at AL = 



27rn 
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L 1 + n'^ir'^a/L 
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i4 



(9) 
(10) 



These analytic approximations give very good agreement 
to the exact numerical solutions of Eq. ([t]) provided one is 
quite high up in the spectrum (optical regime), i.e. when 
n 3> 1, but are only valid in the immediate vicinity of 
an avoided crossing. In particular, being quadratic they 
do not display the asymptotically linear behavior away 
from the avoided crossing that can be seen in Fig. |3] This 
can be remedied by fitting the wave numbers close to an 
avoided crossing to the functions 
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which are the forms suggested by the Landau-Zener 
hamiltonian ( 24 ) , as will be explained in Section |v] The 
gap between the wave numbers of the even and odd global 
solutions at an avoided crossing is given by 2 A /he, where 
for later convenience we have defined A as an energy. 
The parameters 7 and A can be obtained by matching 
the Taylor expansions of Eqns ( 11 ) and ( 12 ) to Eqns ([9| 



7 = 2A fic 
he 

A 



27r^n^a 



L 1 + n^Ti'^alL 



(13) 
(14) 



Note that this expression for A is valid up to second or- 
der in the quantity e = kL — 27rn, i.e. the deviation of 
the solutions from the perfectly localized case, e is small 
if ak ^ 1, which is the case for reasonably reflective mir- 
rors and optical wave numbers. A comparison between 
the values of k found from an exact numerical solution 
of Eq. ([7| , and the approximate analytic solutions given 
by Eqns ( 11 1 and ( 12 ) is included in Fig. |3] 

Away from the avoided crossings the global modes be- 
come localized, each member of a pair localizing on a 
different side of the common mirror. The larger a is, 
the stronger this localization becomes (and the more the 
avoided crossings close up). The general procedure for 
the adiabatic transfer of a light mode from one side of 
the double cavity to the other is now apparent: if the 
common mirror position is initially set to one side of an 
avoided crossing (the optimal position is roughly halfway 
to the next avoided crossing-see Section IV I then the 



localized and global modes approximately coincide and 
the mode is localized on one side on the common mirror. 
If the mirror is now slowly translated to the equivalent 
point on the other side of the avoided crossing, the sys- 
tem will follow the global curve on which it began. On 
the opposite side of the avoided crossing this global mode 
approximately coincides with the localized mode on the 
other side of the common mirror. This procedure is il- 
lustrated in Figure [4] 

It is notable that the wave number of the degenerate 
localized solutions at AL — coincides exactly with the 
lower branch of the corresponding global solutions at the 
avoided crossing, rather than lying halfway between the 
two global branches. This is because the lower global 
branch corresponds to an odd solution in the double cav- 
ity. Odd solutions have a node at the 5-mirror and hence 
do not see the mirror at all. Connoisseurs of the related 
problem of the double well potential in quantum mechan- 
ics [17] will immediately observe that this is the reverse 
of the situation found there, where the lower branch is 
due to the even global solution. This reversal of the or- 
dering of the solutions with respect to the double well 
problem in quantum mechanics is due to the fact that 
for light propagating in ID there can be no evanescent 
waves (assuming the refractive index is > 1 and there is 
no absorption) and so there can be no tunnelling barri- 
ers, only wells. The double cavity system for classical 
light waves is actually equivalent to a quantum reflection 
problem in quantum mechanics. This point is discussed 
further in Appendix [X] 



6 



-a 

3 

e 

03 

aj 

-a 
o 

B 
-a 
-a 
O 



1.0 
0.5 
0.0 

-0.5 
-1.0 
1.0 

0.5 

0.0 

-0.5 

-1.0 
1.0 

0.5 
0.0 
-0.5 
-1.0 
1.0 

0.5 

0.0 

-0.5 

-1.0 
1.0 

0.5 

0.0 

-0.5 
-1.0 







i 




1 








u 
O 



> 



-3 -1.5 1.5 3 -3 -1.5 1.5 
Position X (Units of jim) 



Figure 4. This series of plots shows how each global wave 
mode swaps from one side of the double cavity to the other 
as the common mirror is moved through an avoided crossing. 
The odd mode Un,odd{x) plotted as a function of position x 
is shown in the left column, and the even mode Un,cvcn{x) 
is shown in the right column. Each row corresponds to a 
different value of AL/L: the top row has AL/L — —4 x 10~^, 
the second row has AL/L = -0.6 x 10"", the middle row 
has AL/L = 0, the fourth row has AL/L = 0.6 x 10"**, and 
the bottom row has AL/L = 4 x 10~". Strictly speaking, 
the modes only have well defined parity at AL = 0, but it is 
convenient to retain the labels "even" and "odd" for AL 7^ 0. 
The parameters used to make this figure were L = 1 x 10~"m, 
a/L — 0.3, and n = 128 so that k ■ 



IV. TRANSFER RATIO 

Once the wave numbers have been determined, the 
field ampUtudes A'l^'^ and i?"'*^ on the two sides of the 
common mirror can be calculated. From the continuity 
condition for the field across the mirror we find that 



-An 
Bn 



sin(fc„L2) sin[/c„(L — AL)/2] 



sin(fc„ii) sin[A;„(L + Ai)/2] 



(15) 



where we have suppressed the even/odd labels. Using 
this ratio one can construct the global field distribution 
inside the cavity, as shown in Figure|4]for a pair of modes. 
The mirror only needs to be moved by approximately one 
quarter of a wavelength to achieve the transfer of light 
from one cavity to the other. 
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Figure 5. (Color online) The relative amplitudes An/Bn and 
Bn/An of the modes as a function of the displacement AL of 
the common mirror. The parameters used to make this plot 
were a/L = 0.3, and n — 128. 



The ratio given by Eq. ( 15 1 is a key quantity because 



it measures the degree to which the field is localized on 
one side or the other. In Figure [5] we plot A/B as a 
function of AL, for the case a/L — 0.3. The maxima 
of the ratio A/B are located at values of the mirror dis- 
placement which are roughly halfway between pairs of 
avoided crossings. We shall denote these values of the 
mirror displacement by AL„^. Accordingly, we shall de- 
note the values of the wave numbers at these points by 
fc„^. These wave numbers obey the equation 



cos(fc„*L) + 



sin(fc„iL) f2 + a/L 



k-n-kL 



a/L 



0, 



(16) 



which can be obtained by finding the stationary points 
of the ratio ( jlS] ) with respect to AL, and eliminating 
the quantity dfc„^/dAL, which appears as a result of the 
differentiation, by finding an expression for dkm^/dAL 
from Eq. 

Equation ( |16[ ) can be used to find an analytic expres- 
sion for kni,. In the regime ak ^ 1, the wave numbers 
knif are located very close to the points halfway between 
the avoided crossings, km, w (4n ± l)7r/(2L), where the 



± sign refers to the even/odd mode. Expanding Eq. ( 16 ) 
about these points we find that to first order 



fc„,L = (4n ± 1)V2 - 



2 + a/L 
^■Kna/L 



(17) 



The mirror displacement AL„i, corresponding to wave 
number kn* is found by solving Eq. ([t]) for AL and eval- 
uating it at kn-k, i.e. 



ALni, = -. — arccos 



2 sin fc„*L 



cos fc„*L 



(18) 
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where to obtain the second hne we expanded the arccos 
term under the assumption that akm, ^ 1, and also 
made the approximation sin(A:„^L) « 1. 

Having found expressions for km, and ALm,, we are 
now in a position to obtain analytic estimates for the 
maximum value of the transfer ratio given in Eq. (15 1. 



This ratio can first be simplified by noting that at a max- 
imum we expect the numerator to be near unity and the 
denominator to be near zero. Since it is the inverse ra- 
tio B/A which is large for AL small and positive, let us 
consider that. Expanding the numerator about the point 
where the phase is equal to an odd integer multiple of 
7r/2 gives 

sin (^^{L + AL„)) - ( + ^Ln) - I - Trn 

(20) 

and expanding the denominator about the nearest zero 
gives 



sm I - ALn] 



1 



fc„(L-AL„)-™ (21) 



so that close to a maximum the amplitude ratio can be 
written 



An 



-1 + 1 {knAL + AL„A - (tt + 2mr)Y 



Upon inserting Eqns (171 and (|19|) into this equation we 



still obtain a rather complicated result. However, if this 
expression is expanded for large values of n and ka it 
simplifies to 



a: 



2mT 



(23) 



This result is compared against the exact numerical result 
in Fig. |6] The agreement is excellent, except when a/L 
becomes extremely small, in which case a/c„^ ~ 2mTa/ L 
is no longer large, and the approximations made in the 
above derivation break down. 

As can be seen in Fig. |6j as a increases (i.e. as the 
reflectivity of the common mirror increases), so too does 
the maximum possible shift in amplitude. However, there 
is a trade-off: higher reflectivity gives smaller gaps at the 
avoided crossings, which requires slower mirror speeds if 
the transfer is to be adiabatic. This is discussed in Sec- 
tion |V] Furthermore, whilst it is important to know the 
maximum value A„/_B„ can take in order to evaluate the 
optimal performance of the double cavity as a light trans- 
fer device, it is also useful to know the ratio near AL — 0. 
The reason for this is that, depending on the parameters, 
a sizeable amplitude ratio can already be reached at very 
small displacements, see Figure |5] In practice, therefore, 
it may not be worthwhile going to the maximum, espe- 
cially if that involves a longer transfer time. This issue 
is discussed further in Section |x] Using Eq. ^ or ( 10 ) 
in the ratio (151 yields the desired result. This permits a 
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Figure 6. (Color online) The absolute value of the maximum 
value the relative amplitude Bn/An can take as a function of 
the common mirror parameter a. The exact numerical result 
is plotted as the solid curve, and the approximate analytical 
result, as given by Eq. (23 1, is plotted as the dashed curve. 



Note that the difference between these two curves is hard to 
discern. We set n=128 to make this plot. 



very quick calculation of amplitude ratios near the equi- 
librium position, which along with the results of Section 
[V|may be used to minimize losses. 



V. TIME-DEPENDENCE: THE 
LANDAU-ZENER PROBLEM 

So far we have found the time-independent solutions 
to the double cavity problem: the global modes are the 
static solutions to Maxwell's wave equation for a fixed 
mirror displacement. Even if the system begins in a sin- 
gle (global) mode, a mirror moving at a finite speed will 
induce transitions and the light will therefore be excited 
into a superposition of modes. This is a complicated 
problem to treat analytically, but if we can map the pas- 
sage of our system through an avoided crossing onto the 
related Landau-Zener (LZ) problem in quantum mechan- 
ics, then we can take advantage of existing solutions. 

The LZ problem considers a quantum system with two 
states, the energy separation of which is varied linearly 
in time so that at i = they are degenerate. The two 
states have a time- independent coupling A. Stated in 
this way, the two states being coupled are known as the 
diabatic states, and the time-dependent hamiltonian in 
the diabatic basis is 



Hlz = 



E(t) A 
A -E{t) 



(24) 



where the diabatic energies ±E(t) vary linearly in time: 
E = m/2. (25) 
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Diagonalizing this hamiltonian at each instant of time 
gives the energies E± of the adiabatic states 



E±{t) = ±v/??2t2/4 + A2 



(26) 



which form an avoided crossing. At the point of closest 
approach {t = 0) the magnitude of the gap between the 
two adiabatic energies is 2A, and at t — ±oo the diabatic 
and adiabatic energies coincide. 

The adiabatic states interchange their character be- 
tween t = — oo and t = +oo, which is the crucial prop- 
erty we need in order to transfer modes between the two 
sides of the double cavity. The exact solution of the LZ 
problem shows that if at t — —oo the two state system 
is prepared in one of the two states, then the probability 
that it has made a transition to the other adiabatic state 
by t = +00 is [35ti37] 



(27) 



Plz = exp(-27rAVM) . 



Of course, the LZ problem assumes that there are only 
two modes present, but due to the exponential suppres- 
sion in the adiabatic regime of transitions as a function 
their energy difference, we need only consider a single 
even and odd pair in the double cavity if we start from a 
singly occupied mode and move slowly through a single 
avoided crossing. 



It is tempting to apply the result (27) directly to the 



time-dependent cavity problem by identifying the global 
modes we found in previous sections as the adiabatic 
states 



E^ 



hcke 
hckn 



(28) 
(29) 



and extracting from them the parameters A and d that 
are needed to calculate the transition probability Plz- 
Note that the appearance of Planck's constant in these 
relations is purely for dimensional purposes (we have not 
yet quantized the electromagnetic field, see Appendix \K\ 
for more discussion). In fact, the parameter A has al- 
ready been found above and is given in Eq. ( |14[ ). Fur- 
thermore, if the common mirror is moved at a constant 
velocity v, so that AL = 2vt, then the key requirement of 
the LZ model that the energy difference between the di- 
abatic states is varied linearly with time will be fulfilled. 
From Eqns (11)-(14| we find that d = ^v^. For large 
n this gives 



-d w Sirhcnv/L'^ 



(30) 



which coincides with the result obtained by assuming 
that the wave numbers are those given by localized solu- 
tions. 

However, this whole approach of directly applying 
the LZ results to the time-independent solutions of the 
Maxwell equations is wrong in principle because the 
Maxwell and Schrodinger equations evolve differently in 
time. Therefore, although the time-independent Maxwell 
and Schrodinger equations are identical (compare Eqns 



(A4) and (A5|), one should not blindly apply the LZ so- 
lution to the Maxwell problem. Rather, in Section VI 



we justify the application of the LZ result to the moving 
mirror problem by deriving an approximate version of 
the Maxwell equation which is first order in time and is 
valid in the adiabatic regime. In order to set this up, we 
need to express our system in the LZ diabatic basis. Al- 
though they are similar in character, the diabatic states 
do not correspond exactly to the localized modes unless 
a — >■ oo (the LZ diabatic energies must cross halfway 
between the two adiabatic curves whereas the localized 
modes do not). There is a precise way to obtain the LZ 
diabatic modes and this "undiagonalization" is presented 
in Appendix [P] 

The reason it is advantageous to work in the diabatic 
basis is worth mentioning. Consider the time-dependent 
solution expressed in the diabatic and adiabatic bases. 
Denoting the diabatic modes as 4>l{x) and 0i^(x), we 
have 



£{x,t) ^ AL{t)(j)L{x,t) 

£{x,t)^B^{t)Ue{x,t)- 



- AR{t)(j)R{x,t) 

Boit)Uo{x,t) . 



(31) 
(32) 



The time-dependence extends not only to the superposi- 
tion coefficients but also to the basis modes themselves 
and this complicates the treatment. As can be seen a 
posteriori from the LZ transition probability (27 1, the 



changes in the adiabatic superposition coefficients B^it) 
and Bo{t) are exponentially small during a transfer in 
the adiabatic regime. Almost all the change is incorpo- 
rated in the adiabatic modes. By contrast, the changes in 
the diabatic coefficients during a transfer in the adiabatic 
regime are of order unity because in this basis the system 
changes its state. As pointed out by Zener himself 
the changes in the diabatic basis in the transition region 
(close to the avoided crossing) are typically negligible in 
comparison to the changes in the coefficients, and so one 
can take (f>Lix) and (I)r{x) to be independent of time. 
This approximation is analyzed in Appendix [Ej Thus, 
by working in the diabatic basis we can focus purely on 
the time-dependence of the coefficients AL{t) and AR{t), 
which is a significant simplification. 



VI. FROM MAXWELL TO SCHRODINGER 

Our aim in this section is to derive from Maxwell's 
wave equation a Schrodinger-like (first order in time) 
equation for the electric field amplitudes. This will allow 
us to put the use of the LZ result in the time-dependent 
double cavity problem on a firm footing. 

As discussed above, it is more convenient to solve the 
LZ problem in the diabatic basis. On the other hand, the 
exact solutions we have already found by solving Eq. ([t]) 
are for the global modes, i.e. for the adiabatic basis. We 
shall therefore switch backwards and forwards between 
the two bases as needed. We begin in the diabatic basis 
and substitute expression (31) for the electric field into 
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the Maxwell wave equation d^E = fioe{x,t)df£ , neglect- 
ing the time-dependence of the mode functions (j){x,t). 
The validity of this approximation is examined in Ap- 
pendix |E] This gives 

ALd^cfiL + Ajidl^R ^f^Qsix, t){AL(l)L + Ar^b) (33) 

where the dots indicate time derivatives. In fact, the 
time dependence of the dielectric function e{x,t) implies 
that there are also corrections to the above quoted form 
of Maxwell's wave equation [15]; these corrections are 
discussed in Appendix [F] 

We tackle the left and right hand sides of Eq. (33) 
separately. To simplify the left hand side we substitute 



The definitions (401 and Hll can be combined with ex- 



pressions (D4) and (D5| for cos(6') and sm{6) in terms of 



the diabatic variables to give 



sm{0) cos{e) {loI - ujI) = 2Awav//i 



(44) 



which is one of the combinations appearing in (|38| and 
(39 1. In a similar fashion we find that the other combi- 



the expansions (D6) and (D7| of the localized modes in 



terms of the global modes and then make use of the fact 
that the global solutions, being eigenmodes of the system, 
satisfy 



dlU, = -iioe{x,t)ujlU, 
dlUo = -^iQe{x,t)iAjlUo 



(34) 
(35) 



nations are given by 

sm'^{e)ul + cos^{9)ujI = {E/h - oj^^f + A'^/h'^ (45) 

COS2(0)CJ2 ^ sin2(0)c^2 ^ ^^^f^ ^ ^^^)2 ^ ^2 ^ f^2 (^g) 

Putting all this together, and combining with the right 
hand side of Eq. ([33]) , we arrive at the following expres- 
sion for the two-mode Maxwell wave equations 

2iA/h)L,,,AL^R + {{E/h - Lo^^ f + A^/h^}ALC^L 
+2(A/fi)t^avAj^0L + {{E/h + Lo^,f + A^lh?}AR^R 
= -(Al^l+Arcpi^ . (47) 



to give 



dl(pL = -/io£(a;, t) (- sin{e)LolUe + cos{e)u:lUo) (36) 

(^t'pR — ~Mo£(a^,i) (cos(6')Wg[/e -I- su\{9)ijj^Uo) ■ (37) coupled equations of motion for the amplitudes 



Multiplying through by e{x)(l)L{x) or e{x)(l)R{x) and in- 
tegrating over space in order to make use of the orthog- 
onality conditions (DlOl and (Dll), we can pick out the 



Now that we have made use of the exact results for the 
global modes we can go back to the diabatic basis by 
using (D8) and (D9| to obtain 



Al 
Ar 



2Aa;av/fi 



+ 0fl,[sin(6l)cos(6l)w^ -sin(6l)cos(6l)w^]| (38) 

dl^R = -^io£{x,t){^R[cos^'{9)ul + sm^{^)Lul] 

+ 0L [sin(e') cos(6l)a;2 _ sjj^j-^,) cos(6l)a;2j | (gg) 

In order to explicitly express the adiabatic transfer 
problem in LZ form we let 



2Auj^^/h \ f Al 

(f +-av)^ + #y' U« 
(48) 



The two- mode equations of motion ( 48 1 are second or- 



der in time. They are dominated by the fast evolution 
generated by the diagonal terms containing Wav which is 
an optical frequency of order 10^^s~^. This frequency is 
far greater than the other frequencies in the problem and 
comes from the stationary solutions. We can remove the 
diagonal terms by defining 



= Ve^ + Ayn + Lu,^ 

= -VE^ + Ayh + uj,, , 



(40) 
(41) 



AL/R{t) = AL/R{t)exp 
where 



/ 13L/R{t')dt' 
la 



(49) 



cf. Eqns (111 and ([T2|. In these equations ±E/h are 



the frequencies of the diabatic modes, as discussed in 



Appendix p Indeed, Eqns (40) and (41) can be viewed 



as a way of defining the diabatic modes, the adiabatic 
(global) modes with frequencies U!e and Wo being already 
known from solving the transcendental equation ([t]). We 
have also introduced Wav which is the average frequency 
of the even and odd modes 



f3L/R{t) = ViE{t)/hTio,.y+Ayh^ . (50) 
Time derivatives of AL^R{t) then give rise to the terms 

Al/r i^L/R - '^iPh/RAL/R - iPl/rAl/r - P1/rAl/r) 



X exp 



i / PL/R{t')dt' 



(51) 



Wa 



(42) 



and takes account of the fact that the avoided crossings 
are not located at frequency E/h ^ 0, but are in fact up 
in the optical spectrum. From Eqns ( 11 ) and ( 12 ) we see 
that 

27rnc A 
Wav = —r- + IT ■ 43 



By virtue of the fact that /3l and Pr are so much greater 
than the other frequencies, we make a 'paraxial approxi- 
mation in time' (or equivalently, a slowly varying ampli- 
tude approximation) and drop the first and third terms 
on the right hand side of Eq. (51 1, and thereby obtain an 



approximation to Maxwell's wave equation that is first 
order in time 



dA 



L/R 



dt 



WavA -r 

— Ar/l exp 



±* / Wl - PR)dt' 





(.52) 
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These equations can be further simpHfied by noting that 

E 1 A2 



It 



(53) 



and 



Pn-h^2Eil^\^^^...) (54) 

and truncating each of these expansions at their first 
term. In this way we arrive at 



d {A 



dr V A 



e-^"^ /2 \ ( Al 
e'''"'/2 )\Ar 



(55) 



where r = At//i and ■& = hd j ^ are both dimension- 
less quantities. These coupled equations are mathemati- 
cally equivalent to the time-dependent Schrodinger equa- 
tion id^/'/dr — H{t)iP with the LZ hamiltonian (24), 
as can be seen by letting — > exp[— i?9T^/4] and 
Aji — )■ Anexpli'dr'^ /A]. They therefore justify, within 
their range of validity, the application of the LZ result 
(271 to the problem of a slowly moving mirror in a dou- 



ble cavity. However, it is important to bear in mind that 
these equations come entirely from classical electrody- 
namics. 

Let us examine the accuracy of the approximations 
made in obtaining the Schrodinger-like equations. Con- 
sider a cavity of length L = 100 /xm and a/L = 0.01. 
For the optical wavelength 780 nm we find n = 128 
and Eq. (^14| gives A/h « 7 X 10"s"^ Furthermore, 
_E = A at the avoided crossing and remains of this or- 
der of magnitude. One therefore finds from expansion 
( 53 1 that the first correction to the result /J^, /r — (jJhm is 



of magnitude 3 x 10 "^ojav and so relatively small. Also, 



from expansion ( 54 ) we sec that the first correction to 
the result Pr — Pl — 2£' is of second order and equals 
5 X 10~^E. The smallness of this latter correction is sig- 
nificant because the combination 13b, — Pl appears in the 
phases which are generally more sensitive to approxima- 
tions than amplitudes. 

A further very important check is to see if the neglect 
of second order derivatives in Eq. (51 1 is consistent with 



our final result (55 1. Differentiating Eq. (55 1, we find that 
the largest term neglected, namely Aj^jb, has a magni- 
tude of order (A//i)^^^/fl, whereas the smallest term 

retained, namely 2^/?^/^^^/^, has a magnitude of or- 
der (A//i)a;av^L/fl and so is a factor of 10"* bigger. To 
check this point more systematically, we have numerically 
solved the first and second order equations and their so- 
lutions are compared in Figures [7] and ^ The initial 
conditions were taken to be = 1 and A^ = at the 
initial time t = — 25 and the equations were integrated 
up to the final time t = 25. Note that in order to solve 
the second order equations we also require a second set 
of boundary conditions, e.g. the values of the first deriva- 
tives of the amplitudes at the initial time. Because the 
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Figure 7. (Color online) Comparison of second-order-in- 
time dynamics, as given by the full Maxwell wave equa- 
tion ( 48 1 , versus first-order- in-time dynamics, as given by 



the Schrodinger-like equation ( 55 1 , as the mirror is swept 
through an avoided crossing. The initial conditions are the 
same in every panel: Al = 1 and Ar = at r = —25, 
where r = At/h is dimensionless time. The first-order-in- 
time dynamics are also the same in every panel and are set 
by 1? = M/A^ = 1, which is the single dimensionless param- 
eter appearing in the Landau-Zener formula (271 and in Eq. 



( 55 1. The second-order-in-time dynamics depend additionally 
on the value of A/ftojav; panel a has A/ftojav = 1/100, panel 
b has A/ftojav = 1/500 and panel c has A/fiwav = 1/1000. 
Four curves are shown in each panel: the two curves (black 
and red) that begin on the upper left at r = —25 both give 
Al- The black curve (which climbs highe r) is calculated from 
Eq. (48 1. The red curve is given by Eq. (551. The other two 



curves (blue and green) that end on the upper right at r = 50 
give Ar. The blue curve (which is higher) is calculated from 
Eq. (481 and the green curve is from Eq. (55 1. 
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Figure 8. (Color online) A zoom-in of Fig. [t] combining the 
upper right part of all the panels in order to show the time- 
dependence of An at the end of the sweep. Like Fig. [7] this fig- 
ure illustrates how the solutions of the (second order) Maxwell 
equation reduce to those of the (first order) Schrodinger equa- 
tion as A/(ficJav) — >■ 0. The solid black curve is calculated 
from the first order equations ( |55[), and the remaining curves 
from the second order equations (|48P; A/(/icjav) ~ 1/100 (blue 
short-dashed), A/(fttJav) = 1/500 (green medium-dashed), 
and A/(fia;av) = 1/1000 (red long-dashed). 



values of these derivatives are required only at the initial 
time, i.e. before the mirror has started moving, we can 
take their values from the known solutions for a static 
mirror. We find that the first derivatives of the ampli- 
tudes for a static mirror are given by 



Al 



A 

h 



A 

h 



Al 
Ar 



(56) 



This equation should not be confused with Eq. (55). 



The first order dynamics generated by the Schrodinger- 
like Eq. ( 55 ) depend only upon the value of the combi- 
nation 'd. They do not depend upon the value of A by 
itself. B y co ntrast, the second order dynamics generated 
by Eq. ( 48 1 depend upon the separate values of both d 
and A (or more precisely A//icjav) as can be seen if the 



time variable in Eq. (48) is scaled by h/A, like in Eq. 
(55). Therefore, we expect that as A/ZiWav 0, but 



is held constant, the second order dynamics will reduce 
to the first order dynamics. This is exactly what hap- 
pens, as illustrated by the successive panels of Figure [7] 
In these panels the first order dynamics is held constant 
by fixing 1} = 1, but A/fiwav takes the successive values 
1/100,1/500 and 1/1000. 

The reduction of the second order Maxwell wave equa- 
tions to a first order Schrodinger-like equation has some 
broad implications. In the Maxwell theory the ampli- 
tudes Al and Ar give the electromagnetic energy in the 
double cavity 



1 
2 

£0 



e(x)|f (x,t)pdi 
{\AL\' + \An\') . 



So 



(|i?eP + |SoP) 



(57) 



This energy will not in general be conserved under time- 
dependent boundary conditions. In quantum mechanics 
Al and Ar have a radically different interpretation in 
terms of probability amplitudes and do not by themselves 
give the energy. More discussion of this point is given in 
Appendix [A] Of course, the energy will also not in gen- 
eral be conserved in time-dependent quantum mechan- 
ics, however the hamiltonian structure of the Schrodinger 
equation means that the quantity \Al\'^ + |^_rP = 1 is 
conserved even for time-dependent problems. This im- 
plies that by reducing the Maxwell wave equation to first 
order in time we artificially enforce energy conservation 
in our system. 

Interestingly, whatever the value of A/ZiWav, the agree- 
ment between first and second order in time dynamics al- 
ways seems to be very good at the end of the evolution, 
as can be seen in Figure [8] This is probably connected 
with the quite precise vanishing of the energy difference 
between the initial and final states if they are symmet- 
ric about the avoided crossing, as illustrated in Figure 
[9] This figure shows how energy is initially pumped into 
the electromagnetic field from the moving mirror as the 
avoided crossing is approached and is then removed af- 
ter it has been passed. Indeed, Figure [9] also helps us 
interpret the behaviour of the amplitudes seen in Figure 
[s] whereby and as calculated from the sec- 

ond order equations climb above the values calculated 
from the first order equations (which strictly conserve 
\Al\'^ + = !)• This is a good place to remind the 

reader that the curves plotted in Figures [2] and [3] showing 
the avoided crossings of the wave numbers (or equiva- 
lently, frequencies) do not correspond to energies as they 
would in the quantum mechanical case. Indeed, start- 
ing with Al = 1 and Ar = at t = — 25 corresponds 
to starting approximately in the even global mode on 
the left of the avoided crossing, and providing the mirror 
motion is slow, as here, remaining in that mode. The 
wave number of the even global mode is given by the up- 
per curve in Figure [3] which decreases as it approaches 
the avoided crossing and increases once it has passed it. 
This is opposite to the behaviour of the energy shown 
in Figure [9j and underscores the difference between the 
quantum mechanical and classical optical cases. 



VII. QUANTIZATION 

Everything we have discussed in this paper so far has 
concerned classical electrodynamics. We now wish to 
briefly see how the hamiltonian is quantized in our for- 
malism. However, as emphasized in Section |VI[ we cau- 
tion the reader again that this hamiltonian approach only 
strictly applies for the static mirror case because in the 
dynamic mirror case Maxwell's wave equation provides 
corrections to the Landau-Zener hamiltonian. 

We begin by quantizing in the global basis. The even 
and odd modes provide the normal modes of the dou- 
ble cavity and so the hamiltonian can be written in the 
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Figure 9. (Color online) A plot showing the deviation of 
\Al\'^ + \Ar\^ from unity as a function of time, as calculated 



from the full two-mode Maxwell equations ( 48 1 which include 
the second-order time derivatives. As shown in Eq. ( |57[ ), this 
quantity is proportional to the energy stored jn the electro- 
magnetic field. Each curve in this plot has -3 — 1, but the 
upper curve (blue, short dashed) has A/huja_v ~ 1/100, the 
middle curve (green, medium dashed) has A/ficjav ~ 1/500, 
and the lower curve (red, long dashed) has A/Zicjav ~ 1/1000. 



diagonal form 



H = fkUeblbe + htJoblbo 



(58) 



where b^, is the annihilation operator for a photon in the 
even mode and bo is the annihilation operator for a pho- 
ton in the odd mode. In contrast to the classical case, 
the energy now depends upon the frequency of the light. 
Referring to Appendix [D] we see that we can transform 
to the diabatic basis by letting 



be = cos 9 aj^ — sin 6 
bo = sin ail + cos 6 



(59) 
(60) 



where an and a^ are the right cavity and left cavity an- 
nihilation operators, respectively. We therefore find that 



Ldeb^be + ^ob^bo = (<^e cos^ 9 + u!o sin^ 9)a\^aji 
+ (we sin^ 9 + ojo cos^ 9)a\^aL 
+ cos 9 sin 9 {uo — 0Je)a\^aL 
+ cos 9 sin 6'(wo ~ We)a|^ajj . 



(61) 



Using the explicit expressions for cos 9 and sin 9 given by 
Eqns (|D4|) and ( D5 ) , as well as those for uje and uJo given 



by Eqns (140^ and (41 ), we arrive at the second quantized 



hamiltonian in the diabatic basis 



H — {hLOa.y+E)a\aii + {huJs,v — E)a\aL + ^{a\aL+a\aii) 

(62) 

where ±E are the diabatic energies and A the gap ap- 
pearing in the Landau-Zener hamiltonian (24). 



Writing the hamiltonian in this second quantized form 
makes it clear that A is the matrix element of H between 



the left and right diabatic states, i.e. A = {(j) l\H\(1) ji) . 
When A is small, i.e. when the two cavities are weakly 
coupled, the transition probability per unit time between 
the diabatic states should obey Fermi's golden rule 



dP 27r , 2 

d^ - 



(63) 



where p — L/2Trhc is the density of states in either of the 
uncoupled cavities. This is indeed the case, as can be 
seen by first noting that in the weak coupling limit we 
have, from Eq. (14 1, that 



A 



he 2 he 

mra ka L 



(64) 



where to obtain the second equality we have put k = 
2Tm/L, which is exact in the uncoupled limit. Second, we 
observe from Eq. ( B2 ) that in the weak coupling regime 



ka is large, and the transmission probability can be writ- 
ten 



k^a'^ 



A^ 



(27rpA)' 



(65) 



where we have used Eq. ( 64 1 to obtain the second equal 



ity. Finally, using the fact that dP/dt = T/{L/e), where 
L/e is the "return time" of a photon in one of the cav- 
ities, i.e. the time between collisions with the common 
mirror, we recover Fermi's golden rule Eq. (63). 



We saw in Section |Vl| that in the weak coupling regime, 
i.e. when A/ZiWav — >■ 0, the Maxwell wave equation re- 
duces to a Schrodinger-like equation for the classical elec- 
tromagnetic field. The above quantization procedure 
shows us that in this regime the Schrodinger-like equa- 
tion can be promoted to a true Schrodinger equation. 



VIII. REGIMES OF VALIDITY 

According to the Landau-Zener theory, the adiabatic 
regime occurs when, from Eq. (27 1, 



2^ 

hd 



> 1. 



(66) 



Various physical interpretations can be given to this con- 
dition, as we now show. The rate of change of the sepa- 
ration of the diabatic energies, given in Eq. (30 1, can be 
written 



L/e 



. huj„ 



L/C 



(67) 



where ^Wdop = 2fcti « luj^^v/e is the change in angular 
frequency of a localized mode during the photon return 
time i/c, or, equivalently, the Doppler shift acquired 
when a light wave reflects from a mirro r m oving at speed 
V. Then, using Fermi's golden rule Eq. (63 1 to replace A^ 



with dP/di, we find that the adiabatic condition becomes 



2^ 
M 



dP 



(5wdop dt 



> 1. 



(68) 
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Figure 10. A plot showing the various regimes discussed in 
Section IVIIII as a function of the transmission T of the com- 
mon mirror and ratio of the free spectral range cjfsr to the 
average optical frequency oJav. The Schrodinger-like equation 
( 55 I is valid below the dashed line, which is given by Eq. ( 70 1, 



the left hand side of which we have assumed to be equal to 
10*. This is a very conservative estimate, but is still a factor 
of 10 less than the value given by assuming the parameters: 
L — 100 ^m, a = 0.3L, and A — 780 nm. If the Schrodinger- 
like theory holds, then we can apply the Landau-Zener crite- 
rion for adiabaticity, given by Eq. ( |69| l, and which is plotted 
as the solid line. To plot this line we assumed adiabaticity 
sets in when the left hand side is greater than 10 (the expo- 
nential dependence of Plz certainly ensures that this is the 
case). The regime of adiabatic evolution of the electromag- 
netic field is above the solid line. The shaded region shows 
the regime where the Schrodinger-like theory applies and the 
evolution it predicts is adiabatic. 



This says that the change in frequency of the diabatic 
modes during a photon return time must be much smaller 
than the transition rate between the diabatic states 
(roughly speaking, the rate of photon jumps across the 
common mirror). 

The adiabatic condition must also have a classical in- 



terpretation, because Eq. (271 applies to classical electro- 
magnetic waves. To find this we use Eq. ( |65| to replace 
by the transmission probability T to obtain 



2nA^ _ TuJFSR 



4 Wav V 



(69) 



This says that the ratio of the shift in frequency dur- 
ing the photon return time to the free spectral range 
wfsr = 27rc/L (angular frequency difference between 
cavity modes) must be much smaller than the fraction of 
the incident intensity which is transmitted by the com- 
mon mirror. 

In terms of the same physical parameters as introduced 
above, the condition under which the Maxwell wave equa- 
tion reduces to the Schrodinger-like equation can be writ- 



A 



27r 

T WFSR 



> 1. 



(70) 



This is a distinctly different condition from Eq. ( 69 1 for 



adiabaticity, and this means that it is possible to be either 
adiabatic or non-adiabatic, but still be in the Schrodinger 
regime where Eq. (55) applies, see Fig. [lOl 

It is perhaps counter-intuitive that Eq. (|70p for the va- 
lidity of the reduction of the Maxwell to the Schrodinger- 
like theory does not depend on the mirror speed. This is 
because in this paper we have chosen to work in the di- 
abatic basis, which is not the eigenbasis, except far from 
the avoided crossings. When the dynamics of the elec- 
tromagnetic field is seen from the point of view of the 
diabatic basis it is, therefore, not generally stationary, 
even in the case of a stationary mirror. The characteris- 
tic frequency A//i for the intrinsic dynamics in the dia- 
batic basis is set by the coupling between the two cavities, 
and this can dominate the extra time-dependence intro- 
duced by a moving mirror if the mirror is moving slowly 
enough. Therefore, for a slowly moving mirror the limit- 
ing factor that sets the condition ( 70 1 for the validity of 



the Schrodinger-like theory does not involve v. 

A direct dependence on the mirror speed only enters 
the dynamics of diabatic mode amplitudes in the third 
term on the right hand side of Eq. (51), which we have 
so far ignored. Let us estimate the conditions under 
which we can ignore it. We begin by recalling that con- 



dition ( 70 ) comes from comparing the orders of mag- 



nitude of th e fir st and second terms on the right hand 
side of Eq. (51) for A, namely A w 0[(A//i)2^] and 
PA K, 0[(A//i)a;av^]- The magnitude of the third term 



on the right hand side of Eq. (51 ) is 



^A 



(71) 



The condition under which this term can be neglected in 
comparison to the first order time derivative term, is 



0[/3i//3i] 



AL 



5uj, 



T ■ 



dop 



1 VT 
2^ 



> 1 



.(72) 

This condition is obeyed for most practical situations, 
but does require a certain finite value of the transmission 
T and/or a non-relativistic mirror. If we assume a mirror 
with a/L = 0.3 at a wavelength of 780 nm, then accord- 
ing to Eq. (B2) we have T = 7 x 10~^, implying that the 



Schrodinger-like equation (55) holds if t;/c <C 0.008 



IX. FINITE WIDTH MIRROR 

We now go beyond the (5-function model introduced in 
Section |TT] and consider the case of a mirror with a fi- 
nite width. For many purposes the ^-function model is 
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Figure 11. A comparison of the transmission probability T 
for the finite and 5-function mirrors, as given by Eqns ( |B1[ ) 
and (B2|, respectively. The resonances of the finite mirror 
occur at the peaks of its transmission function. For the pa- 
rameters used in this plot, which are the same as those given 
in the caption of Fig. |12| the two curves intersect near a min- 
imum of the finite barrier transmission. This means that the 
two transmission probabilities remain in good agreement for 
a larger range of k than they might otherwise do. 



perfectly sufficient and has the advantage of being mathe- 
matically simple, but there are phenomena which it is not 
capable of describing such as resonances within the mir- 
ror. In this section we compare the wave numbers of the 
normal modes produced with the 5-function model with 
those of a model of a mirror as a uniform slab of dielec- 
tric material of finite width (we have in mind something 
like the SiN membrane used in [16|, although perhaps 
somewhat thicker). Our main purpose is to check how 
well the (5-function model reproduces the predictions of 
the finite width model. Of course, the finite width model 
considered here does not accurately capture the structure 
of an actual piece of dielectric material, but it does take 
the next step in that direction. 

As a model for a mirror of finite width we take the 
same basic model as considered in Section [ll] except that 
we replace the J-function at a; = by a uniform slab of 
dielectric of width 2M centred at a; = 0. The correspond- 
ing permittivity function is given by 

{eo (1 + h) -M <x< M 
So -Ll < X < -M and M < x < L2 

00 elsewhere 



where = Vl + h is the refractive index inside the slab. 
Imposing the boundary conditions Un{—Li) = UniL^) = 
0, the normal modes can be written as: 

' Asm[k{x + Ll)] ~Li<x<-M 
Unix) = { Bsin[nrkx + (f)i] -M <x < M 
C sm[k{x - L2)] M<x<L2. 



The parallel component of the electric field is continuous 
at the mirror surface, so that provides us with two more 
boundary conditions upon the mode function, namely 
continuity of C/„(x) at x = ±M. Another two conditions 
are provided by the continuity of the parallel component 
of the magnetic field at the mirror surfaces (assuming 
the mirror is non-magnetic so that /i = /iq) which, when 
combined with Faraday's law V x £ = —d'B/dt, lead to 
the continuity of the first derivatives of Un{x) at x = 
±M. Implementing these boundary conditions we find 
that the allowed wave numbers k satisfy the following 
equation 



-n^ sin{2nrkAI) sin 
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(73) 
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= 0, 



which is the equivalent of Eq. ([7| , but for the finite mirror 
case. The mirror resonances occur when 2Mknr — In, 
where I = 1,2,3,..., and so the distance in fc-space be- 
tween mirror resonances is n / {2Mnr). This is much 
larger than the distance between cavity resonances 7r/i 
providing L ^ 2Mnr. 

The transmission T — \t\^ 

by ~~ 



ror, given 



Eq. 



1 



4nJ 



sm^{2nrkM) 



(|B1| 
r 



IS 



of the finite mir- 
Appendix |b] as 

shown as a func- 



tion of wave number by the dotted line in Fig. [TT] These 
are the usual transmission fringes of a Fabry-Perot 
resonator. By contrast, the power transmission of the 
(5- function mirror, given by Eq. (B2) as (l 



a 



2fcV4) 



is the monotonic function of wave number illustrated 
by the solid line in Fig. [TT] The two have the same 
transmission when 



7 2 2 

k a = 



^ ' sin2 {2nrkM) . (74) 

m 2 



In general, this equality is satisfied only over narrow wave 
number regions where the two curves cross, but in Fig. 
11 we have chosen a value of a such that they have simi- 
lar transmissions over a wider range of wave numbers in 
the vicinity oi k = 800. In this region, therefore, the 5- 
function mirror can be a simple model for the more com- 
plicated finite mirror case. Of course, the phase shifts 
produced by the two mirrors are not equal, as one can 
see from Eqns (Bl I and (B2 1, and indeed the finite mirror 



phase changes rapidly in the vicinity of the Fabry-Perot 
transmission peaks. At the transmission minima, how- 
ever, the rate of change of the phase shift with k is much 
slower. Therefore, over the range of wave numbers near 
k = 800 where they have roughly equal values of T, the 



15 




-0.01 -0.005 0.005 0.01 0.015 0.02 0.025 0.03 



Length Difference AL (Units of L) 




-0.01 -0.005 0.005 0.01 0.015 0.02 0.025 0.03 



Length Difference AL (Units of L) 



Figure 12. (Color online) Comparison between the wave num- 
ber structure produced by a mirror of width 2M (solid blue 
curves) and the (5-mirror (dashed red curves) inside a dou- 
ble cavity. This plot is for the case when we are far from 
a mirror resonance. In order to make the comparison we 
matched the transmission probabilities T for each mirror over 
the relevant range of wave numbers as best we could (see Fig. 
111. Up to a relative vertical shift, the agreement is very 
good, including the magnitude of the gap 2A at the avoided 
crossings. In the calculation for the 5-function mirror we set 
a/{L — 2M) = 0.009735. In the finite mirror calculation we 
set M/{L - 2M) = 0.001334 and 2nlM/{L - 2M) = 0.1707. 
Note that to make this plot clearer we have only included the 
wave numbers for a single pair of modes in each case. 



Figure 13. The wave number structure produced by a mir- 
ror of width 2M inside a double cavity. This plot includes 
wave numbers for which there is a mirror resonance. At the 
bottom we see the familiar avoided crossing structure which 
couples even and odd modes with the same index n, but for 
higher wave numbers we pass into resonance and the structure 
changes to one where there are vertical avoided crossings that 
allow for the possibility of adiabatically coupling localized 
cavity modes with different values of the index m. The pa- 
rameters used to make this plot were M/(L-2M)=0.0005005, 
and 2nlM/{L - 2M) = 0.01001. 



EXPERIMENTAL FEASIBILITY 



two mirrors produce phase shifts that do not change ap- 
preciably and are equivalent to small fixed displacements 
of the mirror positions. 

Fig. 12 shows the wave number curves for these two 



cases. Since we are comparing the transmission function 
of a (5-mirror to that of a mirror of width 2M, we have re- 
duced the total cavity length for the (5-mirror case so that 
L — > i — 2A/. This means that the optical path length 
outside the mirror is the same in both cases. Neverthe- 
less, the different phase shifts imposed upon the light 
when it traverses the two different types of mirror pro- 
duces a small vertical shift in the wave numbers. The 
two sets of curves are otherwise almost identical, includ- 
ing the gap 2A at the avoided crossings. 

A very different wave number structure occurs in the 
resonant case when an integer multiple number of half 
wavelengths fit inside the mirror, as shown in Fig. |13[ At 
the bottom of the figure we have the familiar horizontal 
avoided crossings where an even and an odd mode with 
the same value of the index n are coupled. However, as we 
increase the wave number the light approaches resonance 
with the mirror and we find that the amplitude B of the 
light inside the mirror becomes large. Intriguingly, we see 
that the avoided crossings become vertical, meaning that 
we can adiabatically transfer between localized modes 
with different values of n. 



There are many different ways to realize the double 
cavity model discussed here in an experiment. The key 
challenge in all of them is likely to be achieving adia- 
batic transfer before the light leaks out of the system, 
e.g. through the end mirrors, which hitherto we assumed 
to be perfect, but which in reality never can be. Remark- 
ably, this challenge has already been overcome in the ex- 
periment |15| performed in a 6/im diameter ring resonator 
made out of silicon where ultrafast electro-optical control 
of the refractive index allows the effective optical length 
of the cavity to be changed on time scales shorter than 10 
ps, which is the lifetime of a photon in their resonator. In 
this way, the wavelength of light injected into the cavity 
was adiabatically shrunk by up to 2.5nm before leaking 
out. The change in refractive index was achieved by shin- 
ing a laser onto the silicon thereby creating free charge 
carriers which, via plasma dispersion, strongly modified 
its optical properties. The same electro-optical effect can 
also be generated electrically [56] , 

In order to make a double cavity out of this system 
we imagine coupling two ring cavities together. It is not 
necessary to displace a mirror: the difference in cavity 
lengths AL — Li — L2 can be controlled by changing 
the effective optical lengths of each cavity in a concerted 
fashion via the electro-optical effect; see Appendix [C] 
The degree to which the light can be localized on one 
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Figure 14. A zoomed-in plot of the relative amplitudes versus 
mirror displacement AL. It is clear that one only needs small 
displacements to achieve large amplitude transfers. Here 
a/L = 0.3, and n = 128. 



side of the double cavity or the other is set by the di- 



mensionless parameter a/L, as discussed in Section IV 
For a central "mirror" such that a/L = 0.01, the maxi- 
mum amplitude ratio that can be achieved when n = 128 
is approximately |^/-B| ~ 10, and so the probability of 
finding a photon on one side versus the other is 100 : 1. 
For a double cavity of length L = lOOyitm we would then 
require a = lfj,m, and at a wavelength of 780nm Eq. ([S]) 
says that this corresponds to a central "mirror" with a 
reflectivity of i? = 0.94. This degree of coupling is readily 
achieved. 

There are some advantages to using a more reflective 
common mirror. Not only does the maximum localiza- 
tion become larger, but also a significant degree of local- 
ization can be achieved for a much smaller value of AL. 
This is illustrated in Figure 14 for a/L = 0.3 where we 



see that an amplitude ratio of 20:1 can be achieved with 
a displacement of only AL = 10~^L (the maximum am- 
plitude ratio that can be achieved with a/L — 0.3 and 
72 = 128 is 240:1). 

Another possible way to realize the double cavity 
model is to displace the mirrors, e.g. by piezo-electric 
motors. These are routinely used to control mirror posi- 
tions in order to actively stabilize the lengths of ultrahigh 
finesse optical cavities against vibrations [49]. However, 
it is an obstacle that only rather low mirror speeds can be 
achieved. Piezoelectric motors can accurately move small 
distances (nm), but the current fastest ones of which we 
are aware only reach speeds oi v = 0.8m/s [50J. This 
speed can be effectively doubled if all three mirrors are 
moved at once as shown in Fig. [15] On the other hand, 
the distances the mirrors need to move are only a fraction 
of a wavelength: in order to transfer light from one side 
of a double cavity to the other we only need the mirror to 
move by a distance of up to A/4, i.e. AL needs to change 



by A/2. This gives us a transfer time of i = A/4i) 120ns 
if we take v = 1.6 m/s. 

For the case when a/L = 0.01 with a cavity of length 
L = 100/im, the cavity mode corresponding most closely 
to light with a wavelength of 780 nm has n — 128. Insert- 
ing these numbers into Eq. ([l4| gives A//iWav = 3 x lO^''. 
Combining this with Eq. ( ]30[ ), we see that at a velocity 
oi V = 1.6m/s we have M^/A'^ = 2 x 10"** and thus 
the transfer is perfectly adiabatic, i.e. the Landau-Zener 



transition probability Plz given by Eq. (27 1 is utterly 



negligible. Non-adiabatic transitions only become signif- 
icant in the ultrafast electro-optical case: if we assume 
the same parameters as for the Fabry-Perot cavity, ex- 
cept that the transfer takes 10 ps, and thus a rate of 
change of cavity length of dAL/dt — 39, 000 m/s, wc find 
hdn/A^ = 6.8, and a non-adiabatic transfer probability 
of exp [-2t:A'^/{M)] = 0.4. In Appendix |f] we make the 
observation that such a rapidly changing refractive index 
probably necessitates the use of a modified Maxwell wave 
equation for its proper theoretical description. 

At the relatively slow transfer speeds of piezo-electric 
motors the leakage of light out of the end mirrors be- 
comes important. To estimate this effect we note that a 
number of recent experiments have realized Fabry-Perot 
microcavities with lengths below 200/im and finesses of 
up to J" = 10^ [5TI - f53] . The decay rate of the electric field 
amplitude in a cavity is given by k = cn /2dJ- . Being op- 
timistic and assuming that for cavity of length 100/im we 
can achieve a finesse = 10^, one finds k = 5 x 10^s~^, 
and hence the field decays to only exp[— Ki] — 0.6 of its 
original value during a transfer lasting 120ns. Thinking 
in terms of transferring a single photon, the photon has 
an sa 60% probability of escaping the cavity by the end 
the transfer. 

The situation improves when the central mirror is 
made more reflective, as mentioned above. For a/L = 
0.3, which is still an order of magnitude less reflective 
than the end mirrors, and with a total change of the cav- 
ity length difference of AL = 0.0002i (i.e. moving be- 
tween the points AL = iO.OOOlL either side of AL = 0) 
for a cavity of length L = 100/im, the transfer time is 
t — AL/2v = 6ns when v = 1.6 m/s, leading to a 6% 
probability that the photon will escape during the trans- 
fer. 



XI. DISCUSSION AND CONCLUSIONS 

The optical double cavity is a closely related system 
to the double well potential in quantum mechanics, al- 
though, as discussed in Appendix [A| the common central 
mirror in the cavity case does not act as a tunneling bar- 
rier. The wave numbers (frequencies) of the global modes 
inside the double cavity form a net of avoided crossings as 
a function of the difference in length of the two cavities, 
and this permits the adiabatic transfer of light from one 
side of the common mirror to the other. Alternatively, if 
the transfer is halted right at an avoided crossing the light 
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Figure 15. (Color online) By simultaneously moving all 3 mir- 
rors, it is possible to double the effective speed of the common 
mirror, i.e. quadruple the rate of change of the relative cavity 
length AL{t) = Li{t) - L-2{t). 



is left in a superposition of being in both cavities. The 
double cavity is arguably, therefore, the simplest possible 
prototype for a quantum network which can transfer pho- 
tons on demand from one node to another or can place 
them into superpositions of being on different nodes. 

Using a 5-function model for the common mirror, we 
have obtained analytic results such as those given in Eqns 
^ and ([To]) for the wave numbers close to an avoided 
crossing. An expression for the maximum possible trans- 



fer ratio of light is given in Eq. ( 23 1 , and depends on 



the transmission properties of the central mirror, as pa- 
rameterized by the dimensionless quantity a/L. Whilst 
the (5-function model gives an adequate description for 
many purposes, we saw that when the central mirror is 
modelled as a dielectric slab of finite thickness there is 
the possibility of mirror resonances which dramatically 
change the structure of the net of avoided crossings for 
the global wave numbers. In the vicinity of a mirror res- 
onance the horizontal avoided crossings become vertical 
avoided crossings, suggesting a way to adiabatically move 
up the ladder of cavity resonances. 

At first sight, the Schrodinger and Maxwell wave equa- 
tions seem to give incompatible descriptions of the time 
evolution of the electromagnetic field because they arc 
first- and second-order in time, respectively. The first- 
order time evolution found in quantum theory is a conse- 
quence of its hamiltonian structure, which is deeply con- 
nected with the conservation of probability. The stan- 
dard method for quantizing the Maxwell wave equation 
d^£ — {n^ / c^)df £ = 0, where nr(x) is the refractive index, 
is to give it hamiltonian structure: a separation of vari- 
ables of the form £{x,t) — i?„(t)J7„(a;) leads to the 
two equations d^Un + n^Un = and i9^i3„ -I- w^-Bn = 0, 
where = c^/c^ is the separation constant. The first 
equation gives the normal (global) modes, and the sec- 
ond is the equation of motion for a simple harmonic os- 



cillator. This can be transformed i nto ham iltonian form 
by introducing the variables q = ^hj^mioiB* -1- B) and 
p = iyjmhu) /2[B* — B) whose time evolution is given by 
Hamilton's equations x = dpH and p = —dxH, which 
are first order in time. The final step is, of course, to use 
the commutator [x,p] = ih. This prescription is straight- 
forward to apply in cases where the refractive index is 
independent of time, but when it depends on both time 
and space the Maxwell wave equation does not separate, 
there are no normal modes, and a simple harmonic os- 
cillator equation for the field amplitudes is harder to ob- 
tain. The time dependent problem can be tackled in 
number of ways, e.g. by working in a time-independent 
basis. However, for slow time dependence it can be eas- 
ier to work in the adiabatically evolving basis Un{x,t), 
given by the instantaneous solutions to Maxwell's wave 
equation at time t, so that electric field can be written 

£{x,t) = Y,n^ri{t)Un{x,t)exp[-i J* U}n{t')dt']. 

In this paper our primary concern was not the quan- 
tization of the electromagnetic field, but rather the ap- 
plication of a mathematical solution, namely, that of the 
Landau-Zener problem, to the transfer of classical elec- 
tromagnetic field modes in a double cavity using a time- 
dependent mirror (or, equivalently, time-dependent re- 
fractive indices, see Appendix [C|) . In order to apply the 
Landau-Zener model we need a first-order-in-time equa- 
tion of motion. Our approach, which is closely related 
to using the adiabatic basis mentioned above, was to use 
a 'paraxial approximation in time' (slowly varying enve- 
lope approximation) to reduce the Maxwell wave equa- 
tion to the mathematical form of the time-dependent 
Schrodinger equation, given by Eq. ( |55[ ). This reduc- 
tion holds in the regime A/fiwav — >■ and assumes that 
the dimensionless parameters ^ = M/A^, = Ati/h 
and Tf = Atf/h are held constant, is the quantity 
which appears in the Landau-Zener formula (27), and 
ti and tf are the initial and final times, respectively. 
i9 and t appear in both the full Maxwell equation and 
the Schrodinger- like equation, but A/ZiWav only appears 
in the full Maxwell equation. The reduction, therefore, 
does not affect the (lower order) Schrodinger dynamics. 
This process may be likened to the reduction of quantum 
mechanics to classical mechanics by letting h 0, whilst 
holding the classical mechanics (i.e. the action) fixed [55j . 
One consequence of the reduction to the Schrodinger- 
like equation (551 is that energy is artificially conserved. 

In Section IVIII we saw that 



as discussed in Section |Vll 

as a by-product of our approach we were able to quan- 
tize the hamiltonian in the regime A/ZiWav — 0, and 
therefore promote the Schrodinger-like equation to a true 
Schrodinger equation for the quantum amplitudes of the 
electromagnetic field in this regime. One of the lessons 
our treatment of the time-dependent problem teaches 
us is that outside of the weak-coupling regime one can- 
not simply take the solutions from the time-independent 
Maxwell's wave equation and plug them into the Landau- 
Zener theory in order to obtain time-dependent proper- 
ties of the electromagnetic field. 
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There is a simple physical interpretation of the fact 
that the reduction from second- to first-order-in-time dy- 
namics holds in the regime A//iWav ^ 0. In this regime 
the gap is small, which corresponds to weak coupling be- 
tween the two cavities, as is clear from Eq. (65), which 
says that the transmission probability T oc when 
T is small. When the transmission is small the time 
evolution of the amplitudes AL{t) and Aii{t) becomes 
slow, their higher time derivatives can be neglected, and 
hence the reduction to first-order-in-time is valid. As 
pointed out in Section |VIII[ the weak-coupling require- 
ment A / hwa-v — >■ for the validity of the Schrodinger-like 
equation is an intrinsic feature of working in the diabatic 
basis and is not related to the motion of the mirror. The 
weak-coupling requirement could be eliminated by work- 
ing in the adiabatic basis, but then the time dependence 
of the adiabatic basis must be taken into account. 

In this paper we have focussed on the optical regime. 
However, the basic idea of adiabatic transfer of electro- 
magnetic fields between cavities can also be applied at 
other wavelengths, e.g. microwave, which is relevant both 
to micromaser experiments [571 158| , as well as to the new 
field of circuit QED where superconducting qubits play 
the role of "artificial atoms" in microwave resonators 
|59j . Micromaser cavities with Q factors as high as 
Q = 3 X 10"'^'' have been demonstrated fSS" and this gives 
a photon lifetime of 0.2s, relieving some of the problems 
relating to photon decay during the adiabatic transfer. 
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Appendix A: The connection to the double well 
problem in quantum mechanics 

In this appendix we compare and contrast classical 
light waves obeying the Maxwell equations to quantum 
matter waves obeying the Schrodinger equation. These 
differences are well known, but we include them here for 
completeness. In particular, in the main body of this pa- 
per we have alluded to the existence of a close connection 
between the problem of light in a double cavity, and the 
problem of a quantum particle in a double well poten- 
tial. In fact, as we shall see in this appendix, the true 
correspondence is not to a quantum particle tunneling 
between two classically allowed regions through a classi- 
cally forbidden potential harrier^ but rather to a quantum 
particle passing between two classically allowed regions 
separated by a potential well^ all three regions being clas- 
sically allowed. 



Consider on the one hand the Maxwell wave equation 
obeyed by plane light waves with vacuum wave number 
k = Lo/c propagating in one dimension in a medium with 
refractive index nr{x) 



d^ 



+ k^ni.{x)£ ^ 



(Al) 



On the other hand, consider the time-independent 
Schrodinger equation for a particle with energy E = 
h^k"^ /2m and subject to a potential V{x) 



dx"^ 



1 



E 



^ 0. 



(A2) 



An obvious difference between the two is that in the 
Maxwell case the eigenvalue k^ multiplies the spatially 
dependent term ririx), whereas in the Schrodinger equa- 
tion fc^ does not multiply V{x), despite the suggestive 
way we have written it here. Nevertheless, the two equa- 
tions are equivalent if we can identify 



n'^{x) 



1 



V{x)_ 
E 



(A3) 



This equation forms the basis of matter wave optics. As- 
suming that the refractive index obeys > 1, we see 
that the analogy between and light waves and matter 
waves only holds if V{x) < 0. Thus, the potential "bar- 
rier" presented by the mirror in the double cavity is, in 
fact, a potential well if we consider the problem in quan- 
tum mechanical terms. This explains why we find the 
odd adiabatic solution lies below the even adiabatic so- 
lution (see Figure |3|, in contradiction to our experience 
of the double well problem in quantum mechanics. 

In the specific case of the (5-function mirror, which is 
described by the dielectric function e{x) = eQnf.{x) ~ 
£o(l + aS{x)), we have 



d^f 

dx'^ 



+ k^{l + aS{x))£ = 



(A4) 



whilst the equivalent problem in quantum mechanics has 
a particle of energy E subject to a potential V{x) — 
(3S{x), giving rise to the Schrodinger equation 



d^ 
dx^ 



pS{x) 
E 



(A5) 



The two equations can be considered as equivalent if a = 
~f^/E. 

Other differences between the Maxwell and 
Schrodinger cases include the normalization conditions 
and also the interpretation of the wave number as an 
energy. Starting with the normalization, the wave func- 
tion ip in the above equations obeys 1-0(2;) pdx — 1, 
whereas the electric field £ obeys Eq. (3]). This difference 
arises, of course, because is a probability amplitude 
and probability must be conserved, whereas £ has no 
such interpretation. Turning now to the wave number 
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k, if we were solving a quantum mechanical problem, 
the curves in Figures [2] and [3] showing k versus mirror 
displacement AL would become energy curves. No 
such connection exists in the Maxwell case. Indeed, 
in quantum mechanics energy and wave number are 
intimately connected, e.g. in free space E = h'^k'^ /2m, 
but on the other hand a quantum wave has an energy 
independent of its amplitude. The exact converse is 
true for classical light waves: their energy is completely 
independent of wave number, but depends rather on the 
square of their amplitude. A wavelength dependence of 
the energy of light waves only enters the story when the 
light wave is quantized via the photon energy E = hck. 

The problem discussed in this paper is purely linear 
since we have not included atoms in the cavities which 
can couple to the light field, nor have we considered any 
nonlinear optical elements. However, if a nonlinearity 
is introduced into the quantum mechanical double well 
problem then the resulting system has a close connection 
to the Josephson junction |60j . A number of authors have 
exploited this analogy between coupled optical cavities 
(or fibres) with nonlinearity and the Josephson junction 
problem to predict effects such as self trapping of the light 
in the cavities [T^1 - [ST] . If an array of cavities is consid- 
ered, then the possibility exists to realize a bosonic Hub- 
bard model, for which the self-trapping effect develops 
into a full-blown bosonic version of the Mott-Insulator 
quantum phase transition |61[ I62j . 



Appendix B: Transmission amplitudes for lossless 
dielectric mirrors 



In this appendix we collect together some results on 
the transmission functions of (5-function and finite width 
mirrors. The problem of light scattering in one dimen- 
sion from a finite mirror with a homogeneous refractive 
index is essentially the same as that of quantum particles 
scattering from a square potential well |63j (see also Ap- 
pendix |A]). For concreteness, consider an incoming wave 
Aexp(«fca;), a reflected wave i? exp(— ifcx), and a trans- 
mitted wave Cexp(ifca;). The thickness of the mirror is 
taken to be 2M, and it is assumed to be made of a lossless 
dielectric material with refractive index rir- The trans- 
mission amplitude t is found to be 

C exp(— i2Mfc) exp(i arctanf "^''^ tan 2Mknr\) 



Y/l + ^4^sin2(2Mfcn,) 



(Bl) 

The transmission probability T for light is given by T = 
and this is related to the reflection probability R by 
R^l-T. 

Performing a similar calculation for a (5-function mir- 
ror, as specified by the dielectric function e{x) — eo(l + 
a5{x)), we find 



t 



1 



exp(i arctan[fca/2]) 



1 - ika/2 ^H-fc2Q,2/4 



(B2) 



This latter result can also be obtained from the transmis- 
sion amplitude for a finite mirror as given above if we let 
Af — ^ 0, and — > oo, but keep the quantity n^M fixed. 
This procedure allows us to identify the parameter a as 
being given by 



a = 2Mnt 



(B3) 



This only holds providing 2Mknr ^ mn where m = 
0,1,2..., i.e. avoiding resonances whereby an integer 
number of half wavelength fit inside the mirror. 



Appendix C: A double cavity made of two coupled 
waveguides with controllable refractive indices 

In this appendix we consider a double cavity formed 
from two coupled waveguides, each made of a dielectric 
whose refractive index can be controlled, e.g. via the 
electro-optic effect. The optical lengths of the two cavi- 
ties can be independently controlled by modulating the 
refractive index in each waveguide, and this allows the 
optical modes to be transferred between the two waveg- 
uides without moving the mirror. Consider the situation 
where each waveguide has the same length L/2, and one 
has the refractive index ni and the other has refractive 
index n2. The original double cavity model, as depicted 
in Fig. [T] can be converted into the waveguide problem 
by substituting 



^ ^ 2 



«2- 



(Cl) 
(C2) 



To conserve the total optical length during the transfer, 
we choose 



ni 

712 



no 
no 



(C3) 
(C4) 



Here, no is the unperturbed refractive index of the waveg- 
uides and rj is the change due to the modulation. Even 



though the substitutions ( CI I and ( C2 ) capture the main 
relations between the two models, there remain some 
technical differences which originate with the second 
boundary condition given below Eq. We find that 
in the coupled waveguide case the transcendental equa- 
tion Q for the global wave numbers becomes 



tan(n2fcm-^/2) 



"-2 



tan{nikmLl2) 



Til (akm/ni) tSLn{nikmL/2) — 1 



(C5) 

Although the arguments of the tangent functions are ex- 
actly as given by the above substitutions, the factors in- 
volving the refractive indices ni and n2 that appear out- 
side of the tangents are not. However, we note that in 
practice the magnitude of the perturbation rj that is re- 
quired to change the optical length AL — Li — L2 ^ rjL 
by a significant amount can be small. Indeed, according 
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to Section IV complete transfer of a mode is achieved 
when fcAL tt, and if we take A = 780 nm and 
L = 100 /im we then find we require "q = 0.004. This 
should be compared with ng which by definition must be 
greater than unity. To a good approximation one can 
then replace Eq. ( C5 ) by 



tan(n2fcmi/2) 



tan(ni/c„,L/2) 



{akm/nQ)t&n{nikmL/2) - 1 



(C6) 



Note that we have kept the full refractive index depen- 
dence in the arguments of the tangent functions where it 
is crucial. This transcendental equation is obtained from 
Eq. ([6| by precisely the substitutions (CI) and (C2), plus 
the replacement ak — ;> ak/ng. This latter replacement 
can be understood by examining the transmission proper- 
ties of a (^-function mirror that lies between two dielectric 
media, one with refractive index ni, and the other with 
refractive index n2. We find that 



t ■■ 



2ni 



TT-i + ri2 — ika 



so that 



(ni+n2)^ 



1 



4n2 



(C7) 



(C8) 



cf. Eq. (B2|. In the last step we substituted the defini- 
tions (C3) and (C4) and assumed that i] <s:_ uq. Thus, 



when surrounded by a medium of refractive index no, the 
transmission of the i5-function mirror is obtained from 
that in vacuum by putting ak — > ak/no, as above. This 
means that to maintain the same transmission in the two 
cases one should increase either fc or a by a factor of uq. 

With a little care most of formulae given in the main 
body of this paper can be modified for the waveguide 
case. For example, in the key Eqns (11), (12), (13) and 



(14), the speed of light c in needs to be replaced by the 



speed c/no in the guide. 



Appendix D: The diabatic basis 

The solution of the transcendental equation ([t]) allows 
us to construct global modes which exist throughout the 
double cavity. The global modes are the eigenvectors 
of the time-independent Maxwell wave equation in the 
double cavity and need to be recalculated for each value 
of the difference between the cavity lengths AL. These 
modes give the so-called adiabatic basis as specified by 
Eqns (28) and (29), and we denote this basis by Ue{x) 
and Uo{x). 

However, as discussed in Section [Vj it is convenient to 
be able to express the double cavity problem in a ba- 
sis which we refer to as the "LZ diabatic basis", or the 
"diabatic basis" for short. We define the diabatic basis 
as the basis in which solutions to the time-independent 
Maxwell/Schrodinger wave equation can be expressed in 



the form given by the LZ hamiltonian matrix (24). In 



other words, the diabatic energies ±E(t) and eigenvectors 
are those found from the LZ hamiltonian when A = 0. 
We shall denote the diabatic basis vectors 4>l{x) and 

For a (5-function mirror the diabatic basis does not in 
general coincide with the entirely localized basis, as de- 
fined by the wave numbers kn — 21171/ {L ± AL) found 
when a ^ 00. This can easily be seen by comparing the 
LZ hamiltonian ([24]) , for which the diabatic energies cross 
exactly halfway between the two adiabatic solutions, to 
the entirely localized solutions which cross on the lower 
adiabatic curve (see Figure [s]). 

Actually, we are in the slightly curious situation of al- 
ready knowing the adiabatic basis and needing to find 
the diabatic basis. This is the converse of the usual situa- 
tion. Indeed, the adiabatic solutions are in a sense better 
defined because they correspond to the basis which diag- 
onalizes the hamiltonian and is therefore unique. Mean- 
while, the diabatic basis is somewhat arbitrary because 
it corresponds to a particular "undiagonalization" of the 
hamiltonian. Nevertheless, this particular undiagonaliza- 
tion into the specific LZ form ( 24 ) can easily be achieved 



by rotating the states correctly. 

We begin from the diagonalized LZ matrix (i.e. LZ 
hamiltonian in the adiabatic basis) 



ad 



E+{t) 
E^{t) 



(Dl) 



where E± = ±VE^ + A^ are the adiabatic energies 
(global wave numbers), where the plus sign corresponds 
to the solution which is even at AL = 0. We want to 
transform this matrix into the LZ matrix in the diabatic 
basis 



Hlz — 



E{t) A 
A -E{t) 



(D2) 



In order to go between (Dl) and (D2| we perform a 

HlzS, where 



similarity transformation. Let i/ad = S 



S = 



cos 6 sin 9 
— sin 6 cos 6 



(D3) 



We find that the matrix elements of S are given by 



cos 6* — 



and 



sm f = — I 



/ V£;2 + A^ + E 

2V^2 + 



\ 



E< 



El-A^ 



2E. 



(D4) 



' ^E^ + A^-E 



E+ - 



El 



A2 



2E, 



(D5) 

where the expressions after the second equality in both 
equations are written in terms of the adiabatic quantities 
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Figure 16. (Color online) Time derivatives of the diabatic 
amplitudes, as obtained by numerical solution of Eq. (55 1. 
The green curve (that takes the value one at r = t/S./h = —50) 
gives the modulus of dAn/Ar, and the red curve gives the 
modulus of AAl/At. The values of the parameters are set to 
§ = M/A^ = 1. 



obtained by solving the transcendental equation ([t]). In 
particular, the value of A is found from A = -E+|al=o- 

The relationship between the two sets of basis states 
is then 



and conversely 
Uo{x) 



cos 9 Ue{x)+sm9 Uo{x) (D6) 

-sine Ueix)+cose Uo{x) (D7) 

cos^ <f)Fi{x) — sin0 <j)L{x) (D8) 

sin 6* (j)R{x) +cose <j)L{x) . (D9) 



The orthogonality relations satisfied by (j^L/R can be 
found by substituting Eqns ([D8]) and ( [D9| ) into Eq. (|3|. 
One finds 



1 



e{x)(j)L/R{x)(l)L/R{x)dx = 1 (DIO) 
' e{x)(t)L{x)Mx)dx ^ Q . (Dll) 



Appendix E: Time variation of the diabatic modes 

In this Appendix we check the validity of the approxi- 
mation made in Section [Vl] whereby the time-dependence 
of the diabatic mode functions (pL and (j)R is neglected in 
comparison to the time dependence of the diabatic ampli- 
tudes Al and Ar. If the time dependence of the diabatic 
mod es i s taken into account then the right hand side of 
Eq. (33) becomes ALcf)L+'2-AL(t)L+AL(f)L plus an identical 



term wher e L — )> R. The largest terms we are neglecting 
in Section 



VI 



are the cross terms Aj^/rc/ji^/r. 
Even for quite modest values of the reflectivity, the di- 
abatic modes become strongly localized on one side of the 
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Figure 17. (Color online) Early time behavior of the real part 
of Al ■ The sohd red curve gives the numerical solution of Eq. 
( 55 1 and the dashed black curve gives the approximate ana- 
lytic solution given by Eq. ( E6 1 . The values of the parameters 
are the same as those in Fig. |16[ 
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Figure 18. (Color online) Early time behavior of the imag- 
inary part of Al. Two curves are plotted, one giving the 
numerical solution of Eq. ( 55 I and the other giving the ap- 



proximate analytic solution as given by Eq. (E6l, but the dif- 
ference between the two is hardly visible for the times shown. 
Only for times greater than r = —22 do differences really be- 
come apparent. The values of the parameters are the same as 
those in Fig. |16[ 



mirror or the other allowing us to accurately approximate 
the diabatic mode functions by the completely localized 
mode functions given by 



<t>L{x) 



L + AL 



2x 



L + AL 



1 



(El) 



and a similar expression for (I)r. The cross terms referred 
to above can then be shown to lead to new terms in the 
equations of motion (i.e. once the mode functions have 
been integrated out) of the form 



A, 



dx — Atv 



dAL 



(E2) 

and a similar term for (jjR. This needs to be compared 
with Aj^^R. Removing the diagonal terms via the trans- 
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formation (49) to give the slowly rotating amplitudes, we first term on the right hand side) 



find that, assuming AL ^ L, the quantities that must 
be compared are 

-2if3L/BAL/B VS 2{AL/R,-iPL/RAL/R)'^ (E3) 

where the left hand side refers to the existing calculation 
in the main part of this paper and the right hand side 
gives the new terms. In the case when v = 1 m/s and L ~ 
100 /im, the first term on the right hand side is clearly 
negligible in comparison to the other terms. To find the 
magnitudes of the remaining terms let us go back to the 
Schrodinger-like equations ( 55 ) . Taking the derivative of 



these equations with respect to time we find, after some 
manipulation, the following two uncoupled equations for 
each amplitude 



dT2 

d'An 

dT2 



~ —IVT— Al 

dr 

■a difl, T 
= IVT— Ajf 

dr 



(E4) 
(E5) 



where we have defined the dimensionless quantities r = 
tA/h and d = hd / A"^ . In his original paper [35], Zener 
mapped these equations onto the Weber differential equa- 
tion, and used the known asymptotic results for Weber 
functions for r — > ±oo to obtain the transfer probability 
quoted in Eq. (27 1. Here we are interested in a slightly 



different problem, namely the solutions of equations ( E4 ) 
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and (E5) for finite times. In particular, referring to Fig- 
ure 



we see that the time derivative oi A^ takes its 
smallest values at the beginning of the time evolution. 
These are the times we therefore need to worry about 
most since it is then that it is most likely that the mag- 
nitude of the right hand side of Eq. ( E3 1 may exceed the 



magnitude of the left hand side. 

At early times, when Al k, \ and A/j sa 0, we find that 
Equation ( E4 1 is approximately satisifed by 



AL(r) = l + ^ln[T/ro] 



(E6) 



+ X exp[iz9ToV2] {^i[-~i^Tl /2] - Ei[-iz9rV2] 

where Ei(z) is the exponential integral function |64j . and 
Tg is the initial time at which Aj^ — 1. The real and 
imaginary parts of this solution are compared to the nu- 
merical results in Figures (17) and (18). This solution 



yields the following time derivative 



AAl 
dr 



I (l - exp[ii?(ro2 



^)/2] 



(E7) 



In terms of the dimensionless time r, the comparison 
jiven in Eq. ( E3 1 becomes (after dropping the negligible 



dAi 
~d7 



VS Ai 



(E8) 

From Eq. (|E6|) we see that the magnitude of Al is ap- 



'L A 



proximately one, whereas the magnitude of dA^/dr is 
approximately l/ri). Also, in the main part of the text 
we argued that a reasonable value of A//i is lO^^s"^. We 
therefore see that for the case considered in this paper 
(initial times tq = —25 and i9 = 1), the time dependence 
of the diabatic mode functions can be ignored. 



Appendix F: The Maxwell wave equation in a 
moving dielectric 

Leonhardt and Piwnicki |48j have derived a modified 
version of Maxwell's wave equation 



-1 



£(x,t)=0, (Fl) 



in order to describe light propagating in a medium which 
is moving with velocity v. This equation, which is valid in 
dispersionless dielectrics with refractive index nr{x,t) = 
y^e{x, t)/eo, is correct to first order in |v|/c. In our case, 
the moving dielectric medium is the common mirror, and 
our aim here is to estimate the size of the last term on 
the left hand side because we have not included it in our 
treatment. The maximum size it can take is of order 
k'^v/c, as can be seen by taking 2nr(n^ — 1) = 1 and 
letting both derivatives act on the phase (f) = rirkx — 
Lut of the electric waves. To justify our neglect of this 
term we must check that it is much smaller than the 
smallest term we keep in our reduction of Maxwell's wave 
equation to a first-order-in-time equation, as described in 
Section VI Referring to Eq. ( 51 ) , we see that the smallest 
term we keep is /3dA/dt, where /? « Wav In Section VI 
we argue that this term has a magnitude of fidA/At w 



A/h = lO^^s"^ and thus {l3lc^)AA/At 



10 



10 , 



This assertion is supported by the results plotted in Fig. 
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which show that |dA/dT| = (;i/A)|dA/dt| w 1. 

Consider first the case of a common mirror moved by 
piezoelectric motors, so that w « 1 m/s. Then we have 
that k'^v/c « 10^ m^2^ which is considerably smaller than 
the smallest retained term. However, the other case dis- 
cussed in Section [X] concerned a waveguide whose refrac- 
tive index is changed on a timescale of 10 ps, which is 
equivalent to a mirror moving at speeds v « 20, 000 m/s. 
Under these circumstances we expect the modifications 
to Maxwell's wave equation to become important. 
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